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SUMMARY 


An analysis is presented of the load at which yielding first occiirs 
in actual columns, taking adequately into account all the factors which 
have an important effect upon this load. These factors include initial 
defects and the yielding limit of materials. Extensive tests were made 
to verify the assumed relation between the magnitudes of the defects 
and the known properties of columns. The results are expressed as a 
formula or chart applicable to all cases. 


INTRODUCTION 


Investigation of the buckling of columns began in lyklt- with Euler's 
famous theory. Although a large amount of work has been done on this 
problem since that time, the amount of progress from the designer's 
standpoint seems surprisingly small. The classical stability studies 
Initiated by Euler and later extended to cover various types of end 
conditions, variations in cross section, and so forth, consist in the 
determination of the conditions for neutral equilibrium, under infin- 
itesimal displacement, of a perfectly homogeneous elastic column loaded 
along a perfectly straight elastic axis. Classical stability theories 
have been found to be satisfactory for predicting the ultimate strengths 
of "long," that is, very slender, colimms. However, for medivim or short 
columns the defects always present in act\xal columns and the limitations 
to the elastic behavior of actual materials, factors which are not con- 
sidered in the idealized classical stability theories, become of great 
importance. For such columns, which Include most practical applications 
designers still rely upon empirical results expressed in the form of 
curves or formulas, each ciirve or formula being of limited applicability 
These empirical results also determine the range of applicability of 
the classical stability theories and, hence, must be made use of even 
when applying these theories . 
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Buckling problems present certain difficulties by their very nature, 
but the case of the column is the simplest of such problems; and there 
seems to be no very good reason why a rational universal coliamn theory 
should not be developed which would apply equally to all columns and take 
into account all the factors which actually have an important influence 
upon the results. Such a "theory" would, of course, like all theories, 
include a nimber of empirical factors or relations which would have to 
be determined from new or existing experiments; even the classical sta- 
bility theories depend upon the empirically determined stress-strain 
relations of elastic materials. However, the amoimt of en5)irical infor- 
mation required to give such a theory universal applicability would be 
very small compared with what would be required by purely en5)lrlcal 
methods. Such a "universal theory" might be somewhat inconvenient to 
use for design purposes in its complete form, but for the limited ranges 
for which present empirical methods apply it would certainly reduce to 
something of comparable simplicity. The theory could thus replace 
present design methods in these reduced forms even if it were in5)ractical 
for direct iise. 

The advantages of such a development would go far beyond the mere 
replacement of one satisfactory design method by a no more satisfactory 
but more "elegant" method. For example, there is now no way to compare 
one set of empirical results with another set covering a different range. 
Yet, in many fields of engineering such comparisons can be made and prove 
of great value in bringing to li^t and making suitable allowance for 
errors and the effects of variations in testing technique and in the 
interpretations which different investigators put on test results, varia- 
tions which always exist when tests are made and Interpreted by different 
people at different times and places. 

The main advantage of such a development would, however, be the same 
as appears in any field when enpirical results are supplemented by adequate 
general theory. Experimental results are necessarily of limited range. 
Because of the number of variables involved, presently available data on 
colimms - in spite of the great number of tests which have been made - 
cover only a small fraction of possible cases. Only an adequate theory 
can permit safe extrapolation, and the existence of such a theory should 
release designers from design limitations of which they may not even be 
aware . 

Two general criteria are in common use for defining the static 
strength of the parts of machines and structures for design purposes. 

One is based upon the loads at which yielding of the material first 
starts; the other, upon the maximxam loads which can be withstood. The 
first criterion seems logical to use as a basis for design of close 
fitting machine parts which "fail" insofar as serving their purpose is 
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concerned if an appreciable permanent change of shape occirrs. The second 
criterion seems the most logical to use in the design of structures for 
which the exact shape is of relatively little importance compared with 
the ultimate strength. 

Since columns are important elements in both machines and structures 
there should evidently be not one but two column theories, one for the 
colvnnn load at which yielding starts (for which little information exists 
at present) and the other for the ultimate column strength. The present 
paper is intended to supply the first need, namely, a rational analysis, 
supported by tests of a special type, of the load at which yielding first 
occurs in actual columns of any type, taking adeqiiately into accoimt all 
the factors which have an important effect upon this load. 

Although ultimate strengths will not be covered, it is of interest 
here to consider briefly the problem of developing an ultimate-strength 
theory. Up to the load at which yielding starts the action of a column 
is everywhere elastic. Between this load and the ultimate load, part of 
the column is in the elastic state and part in the plastic state (assuming 
that the material has some ductility; if not, the two loads coincide) . 

It is not too difficult to analyze satisfactorily this elastic-plastic 
action for particular cases, and many such analyses have been made; but 
it is much more difficult to set up a general theory covering all columns, 
especially considering the widely varying behavior of different materials 
in the plastic range. 

However, it seems to be general experience that the ultimate strength 
of long columns is only a littxe below the classical stability value, 
while the ultimate strength in the medium range is probably only a little 
above the load at which yielding starts. Only for very short columns, 
approaching something which would usiially be thou^t of as "blocks" rather 
than columns, should the ultimate strength differ very greatly from some 
other known value. Hence, it may be possible to develop a sufficiently 
inclusive ultimate-strength theory by studying in a relatively approxi- 
mate manner the small differences between the ultimate load and other 
known quantities. The dlfficiilty, of coxrrse, is to choose the approxi- 
mations so as to preserve reasonable fidelity over the great range of 
variables required to make such a theory truly "universal." 

This investigation was conducted at the Illinois Institute of 
Technology under the sponsorship and with the financial assistance of 
the National Advisory Committee for Aeronautics. 
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SURW OF PREVIOUS WORK 


Before detailing the present work some discussion should he made 
of previous efforts along these lines (refs. 1 to 9) . While such work 
has shown promise, it has, in the authors' opinion, suffered from 
certain deficiencies which have largely vitiated its \isefulness. The 
distinction between the load at which yielding starts and the ultimate 
load seems to have been given inadequate consideration. Theories have 
been derived for load at which yielding starts and the results of these 
theories have been compared with ultimate load data to determine the 
empirical factors defining the magnitude of expected defects. Where 
direct measurements of defects have been made, they have been confined 
to geometric crookedness j and other kinds of defects, which the present 
tests show to have as great an effect as crookedness, have been neglected. 

The relations which have been assumed between the magnitudes of 
defects and the known properties of the columns also seem both vinreason- 
able and founded upon inadequate dataj it has \isually been assumed that 
defect magnitude is a function of length only or of a cross-sectional 
dimension only or that it is a sum of independent functions of these 
dimensions, whereas certainly the effects of these dimensions are 
actually interdependent and other important factors influence the 
defects. Little thought has been given to putting results in convenient 
general form or to studying such matters as the effects of end conditions 
and variation in cross sections or of the less in^iortant components of 
the defects, all of which must be given adequate consideration before the 
generality of any theory can be considered to be established. 


TESTS 

Specimens 


Because of the large amount of scatter to be expected in the quan- 
tities to be measured - the defects in columns - it was necessary to 
test a large number of specimens. All specimens were tested as colimms, 
and measurements were taken of their deviation from straightness, ini- 
tially and under load. These slender specimens of rectangular cross 
section were made of cold-rolled mild-steel bar stock, cold-rolled 
2024-T5 (24S-T5) aluminum-alloy sheet, and cold-rolled 70T5-T6 ( 75 S-T 6 ) 
aluminimi-alloy sheet; all were of standard manufacture euid cut and were 
handled carefully to avoid introducing any defects not already present. 
Although these specimens were in the long-column range, measurements of 
the second and third harmonics carried the data obtained into the medium- 
column range. 
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To eliminate questions regarding the artificial introduction or 
suppression of eccentricity at the ends, which may arise when hinged-end 
columns are tested - for instance, eccentricities can be introduced 
which add to or partially counteract initial curvatures - all columns 
were tested with built-in ends, as is the case in most practical 
applications. To siinplify the tests and eliminate systematic errors 
due to friction in the measurement of end moments, the tests were made 
with 100-percent end fixity. Such tests, however, are subject to 
systematic errors due to deformations in the specimens or clamps at the 
point of clamping. To eliminate these errors the specimens were held 
in loading heads at some distance outside the points which were taken 
as the ends of the specimens, and small mirrors attached at these points 
detected any rotation, which was then brought to zero by rotation of the 
loading heads. While this system, of course, permits errors, it 
eliminates the systematic errors which might seriously affect the 
statistical information desired. 


Description of Apparatus 

Figure 1 shows a diagrammatic sketch of the loading apparatiis and 
the optical system used for detecting rotations of the ends of the 
effective length of the specimen. The telescope is focused upon the 
image of the scale reflected through the back mirror and small mirror 
on the specimen. With the back mirror placed about 10 feet from the 
apparatus rotations of the small mirror of the order of 0.001° produce 
detectable shifts of the scale point seen against the telescope hairline. 

In the photographs of figure 2, the specimen is shown at (a), with 
the small mirrors defining the effective length at (b) , and with the end 
o lamps in the loading heads at (c) . The load can be measured by the 
dial gage (d) which measures the deflection of the flat springs (e ) } the 
working sections of these springs are machined down from a thicker stock, 
with fillets at the ends, which largely eliminates hysteresis. The 
screw (f) advances the loading head to adjust the axial load, while the 
screw crank (g) rotates the loading head about the axis (h) to bring 
rotations of the small mirrors to zero. 

Deviations from straightness are measured by the micrometer screw (i) 
attached to the carriage (j) which moves upon a track formed of tightly 
stretched piano wires (k) . The micrometer carries a silver-plated tip 
upon its end (l) ) when this tip touches the specimen an electrical circuit 
is completed. By using a galvanometer in this circuit, measurements can 
be made which are accurate to a fraction of a thousandth of an inch. 

During measurement of the deviation from straightness in the initial 
no-load condition, in order to insure freedom from accidental end forces 
and moments, the specimens were held only at the center by a narrow clamp. 


6 


MCA TN 5^15 


Since some of the specimens were very flexible, the weights of the two 
ends of the specimen were balanced by overhead floats at the quarter 
points, as shown in figures 5(a) and 5(b). Nfeasurements of deviation 
were made at the center line of the specimens at eight points along the 
length, as shown in figure 4. 


RESULTS AND DISCUSSION 


As mentioned previously, the inadeq\iacy of classical stability 
analysis lies in the neglection of the limit to elastic action of actual 
materials and the defects always present in actual columns; the defects 
cause bending stresses to develop before the stability limit is reached, 
and these stresses combine with the direct stress (and with any initial 
stress which may be present at the critical point) to precipitate early- 
yielding. 

From the standpoint of column bending the important defects are 
geometric crookedness, lack of elastic homogeneity, and accidental 
eccentricity of loading. All of these have a similar effect in producing 
an initial deviation of the elastic axis of the bar from the strai^t 
line joining the points of application of the resultant axial loads, 
which is called herein the "load line." 

Tn a perfectly homogeneous column the elastic axis, which defines 
the shape of the coliimn for purposes of analysis by classical bending 
theory, would pass through the cen-ters of gravity of cross sections and 
share the geometric crookedness of the outer surface. Because of elastic 
inhomogeneity from slag inclusions, gas bubbles, and so forth, and 
because of the variation in elastic properties in the axial direction 
due to the random orientation of the highly anisotropic crystals of which 
most engineering materials are composed, the true elastic axis will suffer 
an additional deviation from these centers of gravity, passing in effect 
thro\igh the centers of gravity of cross-sectional areas weighted according 
to the local stiffness in the axial direction. Eccentricity of loading 
shifts the load line and thus produces an additional deviation of the 
elastic axis from -this line, as illustrated in an exaggerated manner in 
figure 5* 

For purposes of this investigation all these causes of accidental 
devlationl can be lumped together. This total initial deviation of the 
elastic axis from the load line is designated by the symbol W (as 


^lateral loading and built-in eccentricities also have similar 
effects, and it will be shown that they can be taken into consideration 
along with the defects; however, the latter are the main concern herein, 
since their evaluation is obviously the difficult problem. 
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distinguished from the movement \mder load w) and called herein simply 
the ’‘deviation. ” The starting point of any general column theory must 
he the establishment of laws relating the magnitudes of the important 
constituents of the deviation to the characteristics of columns on which . 
they depend. 

Consider now the best way to measure the deviation and the charac- 
teristics of columns which affect it. 

The deviation W will be some function of the distance x along 
the load line, a different function for each column. The most convenient 
way to describe this function is by the amplitudes of its harmonic 
components, and this proves also to be the most useful way to consider 
its effect upon the buckling process. In the tests, details of which are 
presented in appendix A, the amplitudes of harmonic components of 

the deviation of half wave length were measured over lengths of bar 

corresponding to one wave length of the component. This was done by 
testing lengths of the bars as columns and using an extension of 
Southwell’s method (ref. lO) which had previously been developed in 
reference 11 . A large number of lengths and thicknesses of bars were 
tested; the bars were made of three different standard materials pro- 
cessed by standard methods. As expected, the deviation con5)onents were 
found to depend very much upon the thickness and wave length, the com- 
ponents with larger wave lengths compared with the thickness averaging 
larger in amplitude than the shorter ones. 

Experience has shown that, if a number of similar columns are tested 
which are as nearly identical in every way as it is possible to make them, 
their strengths will vary considerably, but quite definite average and 
limiting (that is, maximum and minimum) strengths can be determined. If 
the deviation components are measured, a corresponding variation (which 
is the chief cause of the variation in strength) will be found, and again 
quite definite average and limiting values can be determined for the 
amplitudes of each harmonic corrqponent. This is what is meant by "average" 
and "limiting" values of such quantities. The variations from the 
average represent true irreducible scatter, which can never be predicted. 
However, the average deviations can be allowed for, and the scatter in 
strength can be allowed for in a more rational and economical way than 
by blanket factors of safety by taking into consideration the maximum 
deviations which produce the minimum strengths. 

If a series of related columns, identical except for a dimension or 
some other characteristic which can be varied continuously, is tested and 
the amplitudes of deviation components are plotted against this 

characteristic, average and limiting ciirves can be determined, which 
describe the function by which the average and limiting values of 

are related to this characteristic. If the relation between the average 
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and limiting values of and all the column characteristics which 

influence them can be determined, proper allowances can be made and 
uncertainty in design can be reduced to true scatter . Insofar as these 
factors are not determined and proper allowances are not made, the 
vincertainty regarding the effect of any characteristic is added to the 
true scatter. 

The characteristics of columns upon which the deviation depends 
might be classified as follows: length and end conditions, size and 

shape of the cross section, the material, and the process by which the 
column is fabricated (which, of course, includes methods of straightening, 
if any, standards of inspection, etc.). The first two, length and end 
conditions, determine the wave lengths which are important in the 
buckling process and, hence, have a very iniportant indirect influence 
upon the deviation^ however, these characteristics are fully taken care 
of if the effect of the wave lengths of the deviation components upon 
their amplitudes is considered. 

The shape of the cross section will usually be associated with the 
fabrication process, and this in txjrn is likely to depend upon the 
material^ these three characteristics are thus closely associated. In 
general, it is impractical to vary these characteristics continuously 
or describe them by numbers. Hence, their effect upon Wjjj, while it 
may be real and important, cannot well be expressed analytically but 
can best be described and taken into account by a nvunerlcal coefficient, 
which is herein designated by C or K and whose value can be tab\alated 
for important distinct combinations of these characteristics. 

Finally, the size of the cross section can, like the wave length, 
be described by a number, and its effect upon Wm can theoretically 
be expressed analytically. For columns of a given shape of cross 
section (that is, for geometrically similar cross sections) the size of 
the cross section can be described equally well by any characteristic 
cross-sectional dimension, such as thickness t, distance from the 
neutral axis to the farthest fiber c, or radius of gyration p (all 
taken for the direction of buckling being investigated) . 

The desired functional relation thus should involve a numerical 
coefficient and three distances Wm, and, say, t. Since it must 

be dimensionally consistent there is no loss in generality if it involves 
only any two independent ratios between these distances, say Wm/t and 
Zm/t. It seems logical to try first a power-f\mction relation between 
these two ratios, which can be expressed as 



( 1 ) 
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where C and n are to be determined. It seems likely that the 
exponent n depends upon broad probability factors and may be substan- 
tially constant for all columns. 

Figures 6(a) to 6(c) show measured values of Wm/t and Im/t 

plotted against each other on a logarithmic scale. Points labeled 1, 2, 
and 5 were obtained, respectively, from the magnitudes of the fundamental 
component and first two harmonics of the total deviation in the test bars. 
Tbp plots show, as is to be expected, a great amount of scatter, but 
they also indicate a definite tendency for Wmy^t to increase rapidly 
as Zm/t increases. The lines marked "max." describe the trend for the 
higher points. The lines marked "av. " should have a somewhat steeper 
slope, corresponding to a larger value of n in equation (l), to fit 
the points best. However, these tests cover the range of wave lengths 
important for medium and long columns but not for short columns. The 
lines shown, when extrapolated into the short-column range, give results 
which are in line with the empirical curves and column formulas in 
common use, while steeper curves would be less conservative^ in the 
absence of data on short columns it seems reasonable to use the relations 
given by the lines shown^ These lines correspond to a value of 2 for the 
exponent n in equation (l) and values of C of about O.OOOO5 for the 
TnavlrmiTTi lines and O.OOOOO7 for the average lines. Even this value of n 
is larger than the values of 0 and 1 which were assumed (on the basis 
of practically no evidence) in the references previously cited, except 
for a recent paper (ref. 12) in which the value of 2 was proposed. 

In the appendix B the following general formula is derived for the 
load upon a column at which yielding starts: 


In this formula Py = ASy is the cross-sectional area A times the 
yield stress Sy (which may be defined in any way desired and reduced 
to allow for initial stresses when this seems justified, as discussed 
in appendix B) , PcZ is the buckling load given by classical stability 
theory (defined as in appendix B in case of a distributed load) , P is 
the correspondingly defined load at which the stress Sy is reached 
at the most highly stressed point, and is the amplitude of that 

harmonic ccmponent of the deviation which has the same half wave 
length Zi as the fundamental (longest) harmonic component of the 
buckling deflection predicted by classical stability theory. The 



(2) 


where 
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length is what has been called the "reduced" or "equivalent hinged 

column" length, so that 


Pc 2 = n2EAp2/zi2 


( 5 ) 


To simplify the final results it is convenient to substitute for 
equation (l) the following equivalent relation: 


cWm ^ 

p2 \ itP / 


_ KZm^ 


2 2 


(n = 2) (4) 


Using this with equations (2) and ( 5 ), the expression for load at which 
yielding starts becomes 


where 


P = N - - PyPci ( 5 ) 



= Py + Pc 2 + KEA (n = 2) 


For some p\arposes it is more convenient to write this equation in terms 
of stresses. Dividing through by the cross-sectional area A gives 


where 


S = 


Q 


- - SySej 


( 6 ) 


2Q= Sy + 



/kE W/2 

\^0ll 


= Sy + Scj + KE 


(n = 2) 


S = p/a is the average stress at which yielding starts, and = Pcz/^ 

rt^Ep2/zi^ is the average stress given by classical stability theory 
(that is, the stress at which instability would occur if the column were 
perfect) . 
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Equation ( 5 ) or (6) can readily te put into a form involving only 
three nondimens ional ratios, say, S/Sy = P/Py, S^j/Sy = arid 

KE/Sy, as follows: 



Sy Sy 


These equations, or other equivalent forms, represent a true 
"universal theory" for column load at which yielding starts. Eqiia- 
tion ( 7 ) can easily be put in chart form; figure 7 shows such a chart 
for the case n = 2, while figure 8 shows how such a chart would be 
affected by different values of n. These charts can be considered to 
be generalizations of the familiar chart of average stress versus 
slenderness ratio and cover the full range from zero to infinite 
slenderness ratio. 

An interesting point brought out by these charts is that only with 
values of n less than 2 would the loads at which very long colimms 
first yield approach the classical stability values. If n = 2 they 
approach values which are equal to Pj,j/j^l + ( KE/Sy . For values of n 

greater than 2 they would approach zero. It is common experience that 
ultimate loads of very long columns do approach the classical stability 
values, but it seems probable from the above that yielding starts at 
considerably lower values. 

Calculations can readily be made from equation (5)j (6), or ( 7 ) or 
charts such as figure 7 , using values of K from tables, of which table I 
may be regarded as a first step; K may be regarded as a "roughness 
factor, " measuring the general roughness of construction. It is a pure 
number, depending upon the associated factors of cross-sectional shape, 
material, and fabrication process; average and limiting values of K 
can eventually be determined for all the combinations of these factors 
of practical importance. This is a large order which, however, it will 
be quite practical to fill in a fairly incl\isive manner by using the 
extensive column data in the literat\jre, that is, by calculating the 
value of K required to make the theory fit such data; these calcvQa- 
tions, however, will have to wait upon the extension of the theory to 
cover ultimate loads, since only ultimate-load data seem to be available. 


12 


NACA TN 5^15 


It is expected that K will not differ widely for variations within 
broad categories such as might be described by the words "refined con- 
struction," "average construction," and "rough construction" and that a 
broad survey of available data, involving the determination of a single 
number to characterize each type of construction, could permit a consol- 
idation of information, with the elimination of many discrepancies, and 
a final relatively simple tabulation from which engineers could choose 
values applying closely to any situation. 

The values of K determined for the small range of column types 
which the present tests cover represent a start in this direction, but 
the main purpose of the tests was to check the general form of equa- 
tion ( 1 ) and determine a reasonable value for the exponent n. As has 
been mentioned, n likely depends upon broad probability laws and is 
subject to little variation. The tests seem to bear out this view. 
Failure to use the most suitable value of n increases the gap between 
the limiting values of K> that is, the proper choice of n is a means 
of reducing unpredictable scatter to the minimum. 


Illinois Institute of Technology, 

Chicago, 111., November 6, 1953- 
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APPENDIX A 

DEVELOPMENT OF HARMONIC ANALYSIS 


From the principles of harmonic analysis, harmonic con5)onents 
Din cos irutx/Z (where m = 1, 2, 3, • • .) S' deviation D(x) oi 
the specimen will have an^jlitudes 


% = 



D(x)cos mjtx/ i dx 


as rfai cos — + bn COS + C-i COS + d-1 

k\^ 8 ^ 8 ^ 8 ^ 


cos + 
8 


9mjt . Ilnur 

d2 cos - — + C2 cos 


1 


In particular, the first three harmonic con5)onents will be 

Dj_ as 0.251 ^a^^ + ag — dj^ — dg^ + O.^l^^bj^ + bg — c^^ — Cg 

D2 ** O.lYY^a^ + ag — bj^ - bg — C]^ — cg + d]_ + dg^ 

D^ as 0. 231 ^.^1^ ^a^_ + ag — d2^ — ^2) ” "^2 ~ ^1 

These formulas permit the determination of the harmonic components 
of the deviation from straightness of the outer surface of the whole 
specimen, initially and under load. In the tests only the symmetrical 
components of the deviation such as those given by equation (9) were 
studied, since the nonsymmetrical buckling modes of a fixed-end column 
are less simple and easy to study by the present methods, and these 
symmetrical components covered as great a range as could have been 
covered by considering the nonsymmetrical modes. For ccmponents such 
as those given in equation (9) it makes no difference whether the 
distances aj^, b^, c^, . . . are measured from the load line or from 
any other parallel or nearly parallel straight line, since a linearly 
varying deviation contains no such components. 



The following definitions are helpful in disciassing the method 
used for determining the total deviation, including the part due to 
inhomogeneities : 
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W’ geometric deviation or crookedness, that is, the initial 

deviation of the median line of the column from the 
load line 

W" nongeometric deviation (due to inhomogeneities); that is, 

the initial deviation of the elastic axis from the 
median line 

W = W’ + W" total initial deviation of elastic axis from load line 


w 


movement due to load 


These are illustrated diagrammatically in figure 9 for the no-load 
and loaded condition of a fixed-end strut such as that used in the tests. 
General expressions for w and W (with similar expressions for W’ 
and W") can he taken as 


w = Wq + Vq x/l + Wjjj cos mjtx/l + > Vp sin prtx/z 

m p 

W = Wq + Vq x/z + "> cos nwx/z + ^ Vp sin pjtx/z 

m p 


( 10 ) 


The mcment equilibrium eqmtion of elementary bending theory is 

-El d^/dx^ = M = Mq + Sq(x + z) + P(W + w) (ll) 


and the boundary conditions are x = +Z 

V = dw/dx = 0 


Substituting expressions (lO) into these equilibrium and boundary condi- 
tions and using the relation = 4rt2El/(2z)^ = n^El/ z^ give 


jwm - [(“i^I'cz/P) - ^'^m| cos mnx/z + 

|vp - pp%cz/P)- ^^p| sin pnx/i + 

Mq/P + Sq(x + Z)/P + (Wo + Wq) + (Vo + Vo)x/z = 0 (12) 



’^o ± ■v’o + - 


Vp sin pjt = 0 



sin nut + pjt 


Vo + 


Vp cos pit = 0 
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These relations are satisfied if, and in general only if, m = 1, 2, 
5 • . p = 1.45^ 2.46, 5-47 • • • and 


Vp - [ 

\ ^ ^1 

0 

) -l]vp=0 

«m - [ 

(m2Pcz/P 

) - ^ " 0 


(15) 


Measurements of a^^, Cq^, . . . and use of equations ( 9 ) in the 
no-load condition give Wjjj', while similar measurements and calculations 
under a load P give Wm' + Wm; subtraction of these gives w^. Knowing 
P and calculating P^j from the dimensions and modulus of elasticity 
of the material, the amplitude of the total deviation components Wju can 
then be obtained from the last equation in equations (I 5 ) . In practice, 
however, it was found easier and more accurate to measure Wm' and two 
values Wjn' + Wn,n and Wj^' + Wj^-^ under two widely different loads Pg, 
and Pj). The term P^^ can then be eliminated between the two relations 

^m “ |^(®^cz/^a) ” iJ^Tna 

% = j^(m2Pg2/Pb) - ij^mb 

giving 

Wj^ = ^nia^mb(i*a ~ ^b) 

'^mai^ “ ^mbi' a 


With this formula for W^, all measurements required are of the same 
type and only relative valioes are needed for the loads Pa and Pjj. 

In figure 10 values of the ratio Wm/t obtained for the 2024-T5 
specimens, are plotted against Zm/t, where Zm = z/m is the half wave 
length of each harmonic component. Points labeled 1, 2, and 5 give, 
respectively, the magnitudes of the fixndamental component and first two 
harmonics of the total deviation, in bars of length 2z. This infor- 
mation is needed in setting up the theory for the buckling of bars of 
length 2Z> the fundamental component is by far the most important 
component, but the higher harmonic components have some effect upon the 
bending stresses produced; and this effect must be evaluated (considered 
in appendix B) before it cein safely be disregarded. 

It was also desired to use the information obtained regarding the 
size of the higher harmonics in order to extend the data regarding the 
size of the fundamental components into the range of shorter columns. 
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This, however, cannot he done directly j that is, the average and limiting 
magnitudes of the third harmonics of hars of length 21 are not neces- 
sarily the same as those of the fundamental components of hars of one- 
third this length. If a har of length 2l is divided into three sections 
and the fundamental con^ionents are determined for each section, then the 
algebraic averages of their three magnitiJdes should he the same as that 
of the third harmonic found from the original har. In many cases, how- 
ever, the fundamental components of the short sections will he of opposite 
sign and will cancel each other as far as the third harmonic of the 
original har is concerned. For the purpose of extending the data 
regarding fundamental components into the range of shorter hars, the 
absolute values of the fundamental components of fractions of the hars 
are needed. These values could he obtained for the geometric devia- 
tion W merely by using known data to make separate harmonic analyses 
for each fraction of the bar. By the same principles as those expressed 
in equation (8) the average of the absolute values of the fundamental 
component of each half or third of a har is 


D2' 


O.lTT^I^'l “ ^1 “ *^i ^1 1 t I a.2 - h2 - C2 + d2 ^ 

0 . 25 ljjo.ll-ll)-(ai - Cl) - (hi - ci)j + 

0.4l4|ci + C2 - d-i - <i2 I + |o.Uli»-(a2 - C2) - (h2 


(15) 


- °2)|] 


The inaccuracy of harmonic analyses based upon so few points is probably 
made up for by the fact that each value of D2' or D5 ' represents an 
average for two or three har lengths. Of coiarse, this averaging process 
also eliminates some scatter, hut the scatter of valiaes obtained from 
such a limited number of points would probably be misleading. 


The second and third harmonics of the geometric deviations of all 
the columns tested were calculated by equations (9) (15) • The 

values for W2' obtained from equation (15) averaged 1.2 times those 
obtained from equation (9), while the values of W^' averaged 2.0 times 
those obtained from equation (9)* Figures 11(a) to U(c) are plotted 
from the data obtained from equation (15); thus, althoii^ the nunibers 1, 
2, and 5 indicate the source of the data, all the points can be taken as 
representing the magnitudes of the fund am ental components of the geo- 
metric deviation of bars. Formxilas (15) could not be used in calculating 
the total deviation, becavise the end conditions of sections of a bar 
under load are obviously not those of fixed ends. However, there is no 
reason to believe that the algebraic averages and the averages of 
absolute values of f\mdamental components of sections of a bar would 
have a different average ratio for bars \mder load (if the sections 
had been tested as separate bars) than for bars under no load. Hence, 
the values for total deviation components W2 were multiplied by 1.2 
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and those for hy 2.0 in plotting the charts of figures 6(a) to 

6(c) and it is considered that these charts therefore show^ to good 
approximation, the magnitudes of fundamental components of the total 
deviation. 

Comparison of the values of the total deviation given hy figures 6(a) 
to 6(c) with values of the geometric deviation given by figures ll(a) to 
ll(c) does not reveal very much difference in slope and in average values 
and not very much difference in the scatter. From this, the important 
conclusion may be drawn that the much easier measurement of geometric 
deviations will hereafter be sufficient and should give results which 
are representative of total deviations. However, this result is in no 
sense due to the effects of inhomogeneities being small - as a matter 
of fact, values of proved to be as large on the average 

as Wjn and Wjjj*, as is indicated by figure 12 for the column made of 

aluminum alloys similar results were obtained with the other two 
materials. The reason why, in spite of this, there is so little differ- 
ence between average and limiting values of W and of W’ is that the 
deviations caused by inhomogeneities W” are as often in the opposite 
direction and subtract from those due to geometric curvature as they are 
in the same direction; hence, these deviations have little effect upon 
the average values and not very much upon limiting values. However, 
it would have been impossible to predict this result in advance or to 
have verified it without an experimental program similar to the present 
one. 
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APPENDIX B 

THEORETICAL DERIVATION 
Simple Case 


Consider first the simplest case, namely, that of a uniform column 
hinged at hoth ends and of such proportions that only buckling in one 
plane, taken as that of the paper, need be considered. Figure I5 shows 
the elastic aocis of such a hinged-end column of length 2 , loaded by an 
axial force P, and with initial total deviation W and movement under 
loeid w. Neglecting the weight or other lateral loading (which can be 
considered by adding the corresponding deflection to W, as discussed 
later) , the equilibrium is given by 

-El d^w/dx^ = P(W + w) 
and the end conditions are 

X = 0,2 

w = d^w/dx^ = 0 

These relations can be satisfied if 

w = ^ Wm sin mnx /2 
m 

W = ^ WjQ sin mjtx /2 
m 

and this expression for W is sufficiently general to represent (that 
is, converge to) any possible deviation shape. Substituting expressions 
(17) into eqiiation (16) and using Pqj = Jt^El/^S gives 

■ [(“^^02/^) " sin nmx/2 = 0 

which is satisfied, in general, only if 

^m = - 1] 


(16) 


(17) 


(18) 
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It is general experience that yielding will first occur in such a 
case because of a combination of the direct stress P/A and the bending 
stress at an extreme fiber of the middle cross section. This will occur 
when the yield-point stress in compression 


Sy = £ + 

A T 


- - 

^ \^Jx=l/2 


= £ - Ec 


P jt^Ec 

z2 

Using relation (18) this becomes 
c = £ + Jt^Ec 

^ A (Pcz/P) - 1 

Multiplying throu^ by A gives 

1 


P = p + Jt^EI ^^1/ 

^ (Pcz/P) - 1 


m^Wjjj sin mjt/2 


m2(Wm/Wi) sin imt/2 
[(n.2Pe,/P) - l]/[(Pcl/P) - 1] 


9W3/W1 


[(9Pcz/!=) - l]/[(Pci/P) - 1 ] 


(19) 


( 20 ) 


Nov, from figure 10 the harmonic components of the deviation of a 
column are on the average related about as 


or 


"fm = cfisr = c(-L)^ 

t \ t / \mt/ 


( 21 ) 


SO that, on the average. 
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Equation (19) then becomes 

cW]_/p^ 


Py - P + Pcj 


(Pcz/P) - 1 


1 + 


(Pci/P) - 1 

5[(9Pcz/P) - 1 ] 


( 22 ) 


The bending stress will increase and approach infinity as P approaches 
Pq2 so that yielding must always occiar before P reaches this value. 

For values of P between 0 and the second term in the braces 

of equation (22) never exceeds For practical struts its maximimi 

value would be considerably less, and further terms of the series would 
be much smaller. Hence, in this case the effect of the higher harmonics 
upon the bending stress can be neglected; in any case they would only 
affect the scatter, since they are as likely as not to cause bending of 
opposite sign from the fundamental component, as suggested by the + sign 
in equation (22) . 

Neglecting all but the first term in the braces of equation (22) 
and solving for P give equation (2), which has previously been discussed. 
It might be pointed out here that in applying relation ( 4 ) to the case of 
hinged struts the values of K found from figure 6 should probably be 
multiplied by 1.2, since, as discussed previously, values higher by this 
amount, on the average, would probably have been obtained haid the funda- 
mental con 5 )onent been measured over half the length (by testing hinged- 
end columns) instead of over the lengths actiially tested. This factor 
has been included in making up table I, so that the values given in this 
table are suitable for hinged-end columns. 


Before finishing with this case seme discussion ml^t be made of 
the effect of initial stresses and lateral loads or built-in eccentric- 
ities. Initial stresses distributed on a microscopic scale (due 
presumably to yielding under previous small loads caused by stress 
concentrations arovind crystals and inclusions) can probably be neglected, 
like these stress concentrations themselves, since such effects are very 
local and scattered and probably have no significant effect on over-all 
shape. However, in cases where significant initial compressive 
stress Si in the axial direction is known to be present on the outer 
fibers of the column (as may sometimes be the case because of rolling 
or other fabrication processes) the stress Si shoiild evidently be 
added to the right-hand side of equation (19) • This is the same thing 
as substituting a "reduced yield stress" Sy' = Sy - Si for Sy, and 
this seems to be the simplest way to allow for such effects. 

Deflections due to known lateral loads and known built-in eccen- 
tricity add to or (just as likely) subtract from those already considered. 
Let the amplitude of the harmonic component of the total deviation due 
to these causes, having the same shape as the fundamental component of 
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the buckling shape predicted by classical stability theory, be 
which corresponds to the amplitude Wp in equations (4) and ( 5 ) due to 
accidental causes. Then, if Kjoax’ is a modified value to be used 
instead, of K ^ny to allow for the effect of Up, there result 


and 





c(Wp + Up) 



Tip 


~ ^Wix 


Kinax' 

give 

J2 . rt^cUpp^' 


7 n 
^1 

it^cUp 


^1^ 



(n= 2) 


(25) 


Thus, should evidently not be changed unless U^_ > W^, in which 

case Kav should be figured from Up instead of Wp. 


General Case 

The foregoing results were derived for the special case of xoniform 
hinged-end columns. It is easy to show that figxares 7 and 8 and the 
equations from which they are derived apply to any column when KE/Sy = 0 
(that is, when there is no deviation of the elastic axis from the load, 
line) provided that or 3^,^ is defined as the classical stability 

limit for the column in question. This is true because when > Py 

yielding evidently will occur as soon as P = Py, while if PcZ ^y 

elastic buckling will occur first but will Immediately result in infinite 
deflections and, hence, infinite bending stresses and yielding, so that 
PcZ = P- 

It is not the purpose of this paper to discuss classical stability 
limits, solutions for which can be foiind in the literature for a great 
variety of columns. The interest here is in the effects of defects and 
of yielding of the material, and it remains to be demonstrated that these 
effects, as exenqplified by the lowering of the curves in flgixres 7 or 8 
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as KE/Sy increases, are the same for all columns. It will he shown 
that with certain simple modifications they prohahly are. Because 
different questions arise in different cases regarding such matters as 
the point in the column where yielding will first occur, it would he 
difficult to set up a general solution covering all types of columns. 

Part of the demonstration will therefore have to he restricted to 
discussion of specific cases; in doing this an attempt will he made to 
span as far as possible the range covered hy actual columns. 

That the foregoing results apply approximately to all types of 
columns can he shown hy the following reasoning. It is well known that, 
whatever the complications - variations in sections, end conditions 
(Including negative fixity), elastic support, and so forth - the 
equllihrium equation for any perfect strut can he satisfied hy an 
infinite number of deflection shapes or "buckling modes," each asso- 
ciated with a particular value of the load. Let w = wif(x) represent 
such a buckling mode where f(x) defines the shape and w^^, the magni- 
tude of the movement, and let Pb he the associated buckling load, that 
is, the load at which equilibrium can exist when the column without 
defects is deflected in this shape. Now compare the equilibrium^ 
equations (representing the equilibrium of external moments and internal 
resisting moments at every section) for this perfect strut and for the 
same strut with an Initial deviation W = Wif(x) having the same shape 
hut a given fixed magnitude defined by 

Then, the term in the equilibrium equation representing the moment 
of the axial force will he P^w = PbWif(x) for the perfect strut, where 
can have any value. For the strut with initial deviation the corre- 
sponding term will he P(W + w) = P(Wy + wp)f(x), where is given 

hut either the load P or w^ is to he determined. All the other terms 
will he Identical in the two equations. Hence, f(x) will also he a 
solution for the second equation (satisfying the same boundary conditions) 
and the following relation must exist between the coefficients of the 
above terms: 


One way to describe the above result is that, if any colu mn has an 
initial deviation in the shape of one of its buckling modes, then a 
movement of this same shape with a magnitude given hy equation (24) will 
occur under an axial load P; this movement tends to infinity if P 
approaches Pb, the buckling load corresponding to this mode. If the 


P(Wi + w^^) = PbW]^ 


Solving for wj^ 



(24) 
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column has an initial deviation consisting entirely of components of 
such shapes^ then corresponding movements given by equation ( 2 ^ 4 -) will 
occur for each of these components and will superpose (assuming that 
the total movements are small) . 


The next question is whether any possible deviation of a column 
can be separated wholly into components having the shapes of the buckling 
modes of the column. It would be easy to answer this question if 
buckling modes were represented by normal functions, like the "normal 
modes" of vibration of an elastic body. Buckling -mode functions are not 
necessarily normal to each other (consider, for example, the symmetrical 
modes of a uniform fixed-end column) , but they nonetheless appear to 
have the property that any possible deviation of a column can be decom- 
posed into components having the shapes of the buckling modes. 


Now, for any end-loaded column (and for any column with loads 
applied between the ends, provided that P and Pb are defined as the 
axial load on the critical cross section due to loads distributed in 
the prescribed manner) , yielding will occur when 


Sy 



ECcr^ll 



cr 





(Pci/P) 



(25) 


where ( d^f /dx^) are, respectively, the area, distance 

to the farthest fiber, and ciirvature at the critical cross section (where 
yielding first occurs), is the lowest of the val\ies of P-jj, and 

f(x) and are, respectively, the corresponding buckling shape and 

the magnitude of the corresponding component of the initial deviation. 
Only this component of the initial deviation is considered in equa- 
tion (25) . The effects of the other components were considered in the 
case of a simple hinged-end column and foimd to be negligible, but this 
must be reconsidered in other specific cases j in any case these effects 
will be negligible when is close to unity, since the primary 

term considered in equation (25) then "blows up" while other terms 
remain small. 


Now, it is general experience that f(x), the shape of a buckling 
mode, is always either a harmonic function (say a sine function but with 
nodes not necessarily at the ends) or close to such a function. If it 
is such a function, equation (25) corresponds exactly to the previous 
derivation and gives the same results. If f(x) is not such a function 
(as in the case of variable cross sections or loads distributed along 
the length) it ceui certainly be represented closely by 
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f(x) = sin rtx/Zi + a sin 2-K-x.fl^ + b sin (26) 

where has the same meaning as in its previous use and where the 

values of a and b are limited by the fact that the c\irvature may 
come to zero (because of large local bending stiffness) in some part of 
the primary wave but cannot reverse in sign. Considering these limita- 
tions, it is possible to calculate limiting values of (d2f/dx2)gj., where 
the critical section is taken as that at which the curvature d^f/dx^, 
and, hence, usually the bending stress, is a maximum. These lifting 
values are from to -2n^/l-^, while the value woiold 

be required to conform to the previous derivation. It can then be con- 

cluded that, if only that component of the initial deviation which has 
the shape of the buckling mode is considered in calculating stresses, 
the results obtained for simple hinged-end colimms can be applied to all 
columns provided that the values of K obtained for simple columns are 
multiplied by factors ranging from 1 to 2. Changes of K of this 

magnitude, of course, produce much smaller changes in P or S, as can 

be seen in figure 7 or 8. 


Extreme Cases 


Now consider some extreme cases more closely . Considering first 
the effect of end conditions, at one extreme there is the case of a 
column free at one end and fixed at the other; this can be considered 
to be half of an equivalent hinged-end column consisting of the column 
and its reflection in the plane normal to the load line at the fixed end. 
If Z is taken as the length of the equivalent hinged-end column the 
entire derivation given previously applies to this case. 


At the other extreme, the case of a fixed-end column has been 
studied previously for a different purpose. The critical section in 
this case will be at one end, where the maximum bending due to the 
deviation components having the shape of the primary (symmetrical) 
buckling mode and the first anti symmetrical mode will add to each other; 
there will always be one side of one end where these and the direct 
stress are all of the same sign, but the stress due to the next symmet- 
rical mode will be as likely to subtract from this, as to add to it. 
Then, using equations (10) and (15), yielding will occur when 


S,, = - 


Ec(^) 

\dx"^/x=0,2Z 


= P + it^Ec 


Wi 


(l/l.E^r(Fc7./F) - 1] . 


(Pcl/P) - 1 - 1] 


(l/^QPclA’) - , . . 

[22(Pc!/P) - J-l 


(27) 
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The maximum value of the quantity in the braces occurs when 

and is 1.5^ ± 0.15. From this it may be concluded that, to be on the 
safe side, the values of Kav obtained for simple colimms should be 
increased by a factor of 1.5^ Q^nd those of 3. factor of IAY- 

However, for Pq^/P — the factor woxlLd be \uiity; for Pqj/P = 2 
these factors would be 1.25 and 1-50, respectively, and so on. 

Consider next the effects of nonuniformity of cross section. Using 
figure 13 let the moment of inertia of cross sections be 


I(x) = I, 


sin nx/Z + ^ a-p sin pjt(x/z) 


sin jtx/Z +^_p^&p sin pjt(x/z) 


(p = 2, 5 , 1^, . . .) (28) 


which^ with siiitahle values of Sp, can readily describe any practical 
variation of stiffness; with one even term 8.2, this expression can 
describe a wide variety of unsyrametrical variations, or with one odd 
term a^, of symmetrical variations of stiffness. Then, the equilibrium 
equation 

-El(x) d^w/dx^ = P(W + w) (29) 

and boundary conditions 

X = 0, Z 

w = d^w/dx^ = 0 

for hinged ends are exactly satisfied if 


w = 

W]_j^sin jtx/z 

+ ^ ap sin pitx/zj 




+ ^ ap sin pjtx/zJ 

(50) 

W = 

W]^ sin rtx/z 


Using expressions (30) 

, equation ( 

'29) is satisfied if 


(jt2/z2)EI< 

jW = P(W + w) 

(51) 

For Wj_ = 0 (perfect 

column) this would become 



*2eIo/!2 = P^, 

(52) 

so from equation (5I) 





Wi = Wi[( 

- 1 ] 
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Yielding will then first occur when 



Ec 


cr 


’£w\ 

^dx^/cr 



rt^Ec 


cr 


Wi 


(Pel/?) - 1 


sin 


jtx/l + 



(55) 


Py - P + PC2 


'cr 


. sin Jtx/l + ^ p^ap sin pnx 


(Pez/P) - 1 


cr 


Since, from equation (52), Iq evidently corresponds to I = Ap^ of a 
unifonn coluinn, it seems very reasonatile to assume that CcrWiA/Io vill 
have about the same average and limiting values as are found for cW^^/pP 
in a uniform column. The equations and charts obtained for uniform 
columns should then apply to nonuniform columns if the values of K 
found for uniform colimins are multiplied by the value of the quantity 
in brackets in equation (55) • 


Choosing x to maximize the expression in the brackets, it is 
found that this quantity may have values as high as 2 for the extreme 
cases contemplated in the previous discussion of the general case. How- 
ever, it is found that this quantity never differs greatly from unity 
for variations which would be used in practical columns and, in fact, 
may be a little less than unity. For columns symmetrical about the 
middle, with a ratio of stiffness in the center to stiffness at the end 
of 2:1, this quantity is about 0.9; for a stiffness ratio of 5 this 
quantity is about unity. The effects of end conditions and of the 
other congjonents of the initial deviation (which are not considered in 
the above discussion) should not be very different from those for uni- 
form columns, and so it may be concluded that the results obtained for 
uniform columns can be applied to practical nonxmlform columns also. 


Next, consider the effect of intermediate loading. For simplicity 
the extreme case was studied of a hinged-end column with axial loads P 
applied at the end x = 0, and at the middle x = z/2 of figure 15. 
Then, a good approximation to the buckling mode should be 

w = W]^jsin rtx/z + a sin 2rtx/ zj (5^a) 

where W]^ and a are to be determined by energy considerations, let 
the component of the deviation of the same shape be 

W = Wp|sin rtx/z + a sin 2rtx/zJ 


(54b) 
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Then, 


the total energy change during a small change in is zero: 



with a similar equation for a small change in a. Using equations (3^a) 
and , carrying out the integrations, and solving the equations 

simultaneously give 


a = 0.066 


1.9 


2^ = P(W + w) 


(56) 


from which 


is obtained. 


Yielding occurs when 




(57) 


The value of the quantity in brackets in eqmtion (57) is about 1.12. 
Thus, it appears that the results obtained for end-loaded coliamns will 
apply closely to columns axially loaded at intermediate points, although 
a small increase in the value of K might be advisable. As mentioned 
previously, P and P^^ must be taken as the axial loads on the critical 
section (where yielding first occurs) due to the load system distributed 
in the prescribed manner. This is necessary in order for the term P/A 
in equation (57) "to represent correctly the direct stress on the critical 
section. In the case just considered the critical section was at 
X « 0.^2 (between the loads) so that this condition was easily satisfied. 
In case of a distributed load the location of the critical section can 
probably be estimated sufficiently closely and should be at about the 
above location for any anti symmetrical distribution of load such as a 
uniformly distributed load. 


As a final example, consider the case of a column on an elastic 
foundation. If the hinged -end column of figure 13 has an elastic 
support of p (force per unit length per unit deflection) the equilib- 
rium and boiindary conditions can be satisfied by the same expressions 


28 


MCA TN 5415 


(eqs. ( 17 )) for w and W as were used for the simple col\mm. In a 
previous paper (ref. 11) it is shown that in such a case 


\ Z2 


£i_\ 




m 




( 58 ) 


vhere the expression in the parentheses on the left side represents the 
buckling stress corresponding to the buckling-mode shape sin nutx/z, 

as can be seen by letting Wj^ = 0 in equation ( 58 ). Buckling will occur 
with the niJinber of half waves m = m' corresponding to the smallest 
value of Pm» that is. 


P 


Cl ~ 





(59) 


The smallest value which the expression in the parentheses on the right 
side can have is 2^ which occurs when 

/il ^ /ll = 1 (40) 

Z2 P h 

for which 

Pcj = 2 /Si 


Assuming for simplicity that the length is such that m' given by 
equation (ko) is an Integer, and considering only three components of 
the initial deviation W, corresponding to m', m’ + 1, and m' - 1, 
equations (l?), ( 38 ), (40), and (2l) lead to the following condition 
for yielding: 


% 




jt^Ec 

(l/m)2 (Pc^/P) - 1 


sin ± 




T(m' + l)2rt2EI 1 i2q 

/pi 

_ I 

z2 (m' + l)^^_ 

/ J 


L - Xp - i)H"' - 

:/z] 

fRm* - l)2jt2EI . z2p ~| 

/pi 

1-1 

,2' '(m'-l)2«2j 

/ : 

1 , 


m 
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For Pcz/P — value of the largest term in parentheses in equa- 
tion ( 4 l) is unity and the results derived for simple struts then apply 
exactly to the present case. The expression has its largest value if 
Pc 2/P — >00 when it becomes 


sin 


m*jtx ^ 
I 


sin 2 sin 

m'+l I ,m'-l Z 


fm' + l\2 


/m' - 1 


V m' / 

Im' + 1/ 

\ m* 

) ^ 1 


( 42 ) 


This is a fixnction of m' and can he evaliiated without great difficulty. 
It is found that points can always he found where stresses due to the 
first and second, or the first and third, terms have the same sign as 
the direct stress" (whichever combination gives the largest value is 
chosen to determine K^v ) f ^ut the stress corresponding to the remaining 
term will then be as litely to subtract from as to add to this amount 
(and so may be used to determine Knja,x) • 

Using these findings it is determined that for 'Pci /^ — "the 
value of Fav for sin5)le columns should be multiplied by approximately 

. - [x/(..) 3 ] and the value of K^ay , by 5 - j^2/(m')2j, where m' is 

the number of half waves in which the column will buckle according to 
classical stability theory. Thus, a column with a mild elastic support, 
such that the buckling shape is still one half wave, would require no 
correction for K. However, a coliamn which is so long that it buckles 
in many waves but is at the same time so strongly supported that P is 
small compared with the classical stability limit (that is, the colijmn 
is in the short-column range) may require corrections by a factor as 
high as 2 for K^y and 5 for K^nY, because there is only a slight 
difference between the classical stability load and the loads corre- 
sponding to neighboring buckling modes; and, hence, the corresponding 
components of the deviation contribute a good deal to the total 
bending stress. 
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TABLE I 

VALUES OF ROUGHNESS FACTOR K 

Values are for simple columns but may be used for other 
cases as discussed in appendix B 



^^D3ax 

Kav 

Standard cold-rolled steel bar stock 

0.00015 

0.00005 

Strips cut from standard flat sheets of 
2024-T3 and 7075-T6 alimilnum alloy 

.00019 

.00004 

Columns of "refined" construction 
(tentative values) 

.00015 

.00005 

Columns of "average" construction 
(tentative values) 

.00040 

.00010 
1 
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( a ) Loading head . 



(b) Deflection-measuring unit. 
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Figure 2.- Photographs of test apparatus. 
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Figure 5*- Sketch showing method of balancing weight of slender specimen. 
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Figure 4 .- Sketch showing points at which initial geometric deviations were measured. 
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Figure 5.- Shift in load line produced by eccentricity of loading. 
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(a) Cold-rolled mild- steel bar stock. 
Figure 6 .- Total initial deviations of columns. 
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(b) 2024-T3 aluminum alloy. 
Figure 6.- Continued. 
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(c) 7075 - t 6 al\aminum alloy. 
Figure 6 .- Concluded. 
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Figure 7*- Universal column chart of 
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Figure 8 
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Figure 9.- Diagram illustrating factors involved in total deviation 
and movement for a fixed-end strut. 
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Figure 10.- Harmonic components of deviation. Specimen of 2024-T5 aliamin\am 

alloy. 
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(a) Cold-rolled mild-steel bar stock. 


Figure 11 .- Geometric deviations of columns. 
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(b) 2024-T5 alTominuin alloy. 
Figure 11.- Continued. 
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(c) 7075- T6 aluminiiiii alloy. 


Figure 11 .- Concluded. 
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Figure 12.- Nongeometric deviation. Specimen of 7075“T6 aluminvim alloy. 
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Figure 15.- Uniform column hinged at hoth ends. 
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